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Abstract 

With an integrated luminosity of 2.47 fb”^ recorded by the ATLAS experiment at the LHC, 
the exclusive decays and ^ Jlij/K*^ of B mesons produced in pp collisions 

nt ^fs = 1 TeV are used to determine the ratio of fragmentation fractions fs!fd- From the 
observed ^ and ^ yields, the quantity ^ measured to 

be 0.199 + 0.004(stat) + 0.008(sys). Using a recent theory prediction for yields 

^ = 0.240 + 0.004(stat) + O.OlO(sys) + 0.017(th). This result is based on a new approach 
that provides a significant improvement of the world average. 
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The production rate of mesons is a product of the bb cross section, the instantaneous luminosity 

and the probability that the l^-quark is bound to an s-{d-) quark. The latter, denoted by the fragmentation 
fraction fs {fd), depends on the probability that in pQCD-inspired calculations [1, 2] a soft gluon splits 
into ss (dd) and that the overlap of the b and s (d) wave functions is sufficiently large to produce a 
bound state. In a similar fashion, B^ mesons. Be mesons and Z>-baryons are produced at the LHC 
with respective fragmentation fractions fu, fc and /baryon- The fragmentation fractions are about 40% each 
for u- and r/-quarks, 10% for 5-quarks, at the percent level for c-quarks and ~ 8% for baryon production 
satisfying the constraint fu + fd + fs + fc + /baryon = 1 ■ Precise knowledge of the fragmentation fractions 
is essential for measuring Zj-hadron cross sections and branching fractions at the LHC. In particular for 
rare decays, such as the branching fraction measurement of ^ [3-5], a precise knowledge of 

fs!fd is important since it improves the sensitivity of searches for new physics processes beyond the 
Standard Model (SM). The fragmentation ratio fs!fd is a universal quantity that was measured by LEP 
experiments [6], CDF [7] and LHCb [8, 9]. This Letter presents a measurement of fs! fd using ^ Jiff 
and ^ JjfK*^ decays. 

The ratio of fragmentation fractions fs! fd is extracted from the measured — > Jjff and ^ JjfK*^ 
signal yields, Mgo and N^o. These are converted into meson yields after dividing by the 

branching fractions of the relevant decays and correcting for the relative efficiency “Reft that is expressed 
as a product of acceptance and selection efficiency ratios for the two modes and is determined from Monte 
Carlo (MC) simulations: 


fs _ ^ ^ K-n-) 

fd ^ J/ff) m ^ ^ ^ 

d 

where the Jjf, f and K*^ are reconstructed in their J/f —> ^^IJ~, f —> K^K~ and K*^ K*7:~ final 

states [10], respectively. The data sample consists of pp collisions collected with the ATLAS detector at 
- V TeV corresponding to an integrated luminosity of 2.47 + 0.04 fb“/ The ATLAS multipurpose 
detector is described in detail in Ref. [11]. 

The PYTHIA 6 and 8 [12, 13] MC generators with parameters tuned to reproduce ATLAS data [14] 
are used to simulate background and signal events, respectively. For the signal channels, the angular 
distributions are produced with the measured polarization parameters [15]. The detector response for the 
generated events is simulated with Geant4 [16, 17]. 

The ^ Jjff and —> JjfK*^ signal candidates consist of two muons and two hadrons originating 
from a common secondary vertex. The Jjf candidates are selected from the dimuon trigger sample re¬ 
quiring two oppositely charged muon candidates, each having a transverse momentum of px > 4 GeV. 
Reconstructed muon candidates are categorized either as combined or segment-tagged muons. A com¬ 
bined muon consists of an inner detector (ID) track combined with a muon spectrometer (MS) track using 
tight matching criteria, while a segment-tagged muon requires an ID track and track segments in the MS 
that are not reconstructed as an MS track [11]. The two muons, of which at least one must be a combined 
muon, are fitted to originate from the same two-track vertex. The vertex fit chi-square per degree of free¬ 
dom (dof) is required to be;t^^/dof < 10. To improve the sample purity, each muon track must have at 
least one hit in the pixel detector, more than five hifs in fhe silicon strip detector and at least one hit in the 
transition radiation tracker that reduces the pseudorapidity coverage to |p| < 2.0 [18]. 

Since the dimuon mass resolution is different for muons reconstructed in the endcaps (1.05 < H < 2.5) 
and for muons reconstructed in the barrel dpi < 1.05), all accepted Jjf candidates are divided into three 


2 



classes: two barrel muons (BB), one endcap and one barrel muon (EB), and two endcap muons (EE). The 
parameters describing the dimuon mass distribution in the 7/i/^ signal region for the three pseudorapidity 
classes in data and in J 1^(1) and —> JI MC signal samples are extracted from maximum- 

likelihood fits. Signal events are selected requiring mass windows of ±3cr around the Jj^ peak in data 
and simulations. Eor data, the selected signal regions are 2.991-3.197 GeV for BB, 2.955-3.235 GeV for 
EB and 2.914-3.275 GeV for EE classes, while in simulations they are slightly smaller. 

The candidates are reconstructed from a T/i/^ candidate plus two oppositely-charged hadrons with a 
kaon mass hypothesis assigned. The dimuon mass is constrained to the Jj^ mass [15] and the //i/^ and 
two kaons have to originate from the same vertex. All combinations are accepted if > 8 GeV, 

l^oi < 3 for the vertex ht and the K^K~ invariant mass lies in the range determined by +2 natural 
widths (T^) around the (p mass peak, 1011 < mx+R- < 1028 MeV. The niR+R- distribution is modeled 
with a Breit-Wigner line shape convolved with a Crystal Ball function [19]. The selected mass window 
retains 85% of signal events. 

The candidates are reconstructed in a similar way. Here, one track of the K*^ decay is assigned 
a kaon mass hypothesis and the other track a pion mass hypothesis. Since ATEAS has limited kaon- 
pion separation capability in the momentum range relevant for this analysis, both Kn mass assignment 
combinations are tested. That with mass closest to the nominal mass is chosen yielding the correct 
Kn selection for 86% of all K*^ candidates. The probability density function (PDE) for the invariant mass 
of correctly selected Kn candidates is modeled with a relativistic Breit-Wigner line shape convolved with 
a Crystal Ball function, while that where the K and n are swapped is modeled with a Gaussian function. 
The decay ^ produces a peaking background in ^ that appears in the low K*^ 

mass region. To remove this contribution, the selected K*^ region is constrained to one K*^ decay width 
around the K*^ mass peak, corresponding to 847 < iurj^ < 942 MeV for data. Since the K*^ line shape is 
narrower in the MC simulations than in data, the Kn mass selection needs to be adjusted in simulations to 
produce identical efficiencies in data and simulations. Eor the K^K~ mass selection, a similar procedure 
is used. 

The signal-to-background ratios for ^ and ^ Jjij/K*^ decays are optimized using three 

variables with high background suppression power: the;^f^/dof of the B vertex ht, the transverse decay 
length Lxy dehned as the length of the vector from the primary vertex (PV) [20] to the B decay vertex 
in the transverse plane, and the pointing angle a dehned as the angle between the B meson transverse 
momentum and L^y If more than one PV candidate exists, the one is selected for which the sum of 
squared transverse momenta of all tracks originating from the vertex, ^ p^, yields the highest value. The 
X^ jdof, Lxy and a selection criteria are optimized using simulated —> Jlipcp and —> Jj^K*^ events 
for signal and data sidebands for background. 

To produce similar and p distributions in data and MC, data-driven weights are obtained by the 
following procedure. Sideband-subtracted ^ Jlt//K*^) pj and p distributions from 

data are compared with corresponding distributions in simulation in the signal region, 5.32 < < 

5.42 (5.21 < ifij/^i^Rta < 5.35) GeV. The upper and lower sidebands 5.20 < mji^^ < 5.25 (5.09 < 
niji^K,Q < 5.16) GeV and 5.48 < < 5.53 (5.40 < mj/^R.o < 5.47) GeV are selected such that their 

summed yields represent the expected backgrounds in the signal region for the data. The weights are 
obtained by dividing the yield in each pj and p bin in data by the corresponding yield of the MC sample 
using only events with odd event numbers. Thus, for each bin (/) and (j) of the pj and p distributions, a 
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weight is determined as a product of a px-dependent and 77 -dependent weights: 


Wij{pT:,r]) 




( 2 ) 


where is the normalized number of entries in the pj bin i and is that in the 77 bin 

j. To obtain good agreement between data and simulation, the procedure is repeated twice. The two sets 
of weights are multiplied and are used to correct the pj and 77 distributions of the MC sample with even 
event numbers. From the corrected MC samples, distributions for;^^/dof, Lxy and a are determined that 
are in good agreement with those measured in the data. The correlation between and 77 is small and is 
accounted for in the systematic error. 

For both modes, the dominant background originates from a 7 / 1/7 produced at the PV plus two oppos¬ 
itely charged hadrons (direct 7 / 1 / 7 ) [21]. Since the hadrons are not associated with any decay, the 

Jlij/K^K~ invariant-mass spectrum does not peak but decreases with mass. Another large 

background consists of two random low-momentum, oppositely charged muons combined with two ran¬ 
dom charged hadrons. Here, the dimuon mass distribution does not peak at the 7 / 1/7 nor does the four- 
particle mass show any peaking structure. Inclusive decays B JfifrX where X is a single hadron or 
a collection of hadrons provide a source of background that is very similar to the signal. If X consists 
of exactly two charged-particle tracks (without any ;r°), the mode is topologically indistinguishable from 
the signal mode. Self-cross-feed in which one or both hadrons from the </> (X'*^) decay are replaced with 
random hadrons is negligible. In addition, peaking backgrounds from ^ and ^ 

contribute to —> Jltf/(p while —> JlipK^n~ also contributes to ^ Jli/zK*^. 

To reduce these backgrounds, the;y^/dof, Lxy and a selections are optimized for each mode separately by 
determining the maximum value of S / + B as a function of selected values for the observable to be 

optimized, where S represents the signal yield obtained from simulation and B is the background extracted 
from data sidebands. For the (B°) mode, the optimization yields;y^/dof < 2.4 (2.6), Lxy > 0.26 (0.30) 
mm and a < 0.14 (0.12) rad. In combination with the 7 / 1/7 mass requirement, the;y^/dof selection reduces 
the combinatorial background significantly, while the Lxy and a selections remove most of the direct 7 / 1/7 
background. 

In the final sample, fhe signal yields Ng) and Ngo are exfracfed from unbinned exfended maximum- 

d 

likelihood fifs fo fhe JliJ/K'''K~ and invarianf-mass specfra, respecfively. The signal PDF is 

modeled wifh fhree Gaussian funcfions wifh common mean fhaf is defermined from fhe fil while widfhs 
and fraclions are fixed fo fhe values obfained from MC simulafions. To accounf for possible widfh dif¬ 
ferences in fhe fwo narrowesf Gaussian funcfions befween dafa and simulafion, an addifional scale faclor 
is infroduced, which is leff free in fhe fif. The peaking background PDF is modeled wifh a Crysfal Ball 
function wifh paramefers fixed fo fhe values obfained in simulafions. The peaking background yield of 
652 ± 93 evenfs is calculated from fhe signal yield. The selection efficiencies of bofh peaking back¬ 
ground modes are defermined from simulafion and are fixed in fhe fif fo dafa. The remaining residual 
backgrounds are modeled wifh an exponenfial function leaving fraction and exponenf free in fhe fif fo 
dafa. 

The signal PDF is paramefrized wifh fhree Gaussian funcfions fhaf describe bofh fhe correcfly re- 
consfrucfed and swapped K^n~ evenfs. The PDF of fhe peaking background is modeled wifh a sum of 
Crysfal Ball and Gaussian functions for which fhe relative —> JliljK^n~ yield wifh respecf fo fhaf of fhe 
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Figure 1: The invariant-mass spectra of (left) and —» Jjij/K*^ decays (right) for data (points with 

error bars), total fit (solid), signal (dashed), residual background (yellow shaded), partially reconstructed events 
(magenta shaded) and peaking background (blue shaded). 

^ JI signal is determined from the corresponding branching fractions and selection efficiencies, 
yielding (4.7 ± 2.4)%. Most of the residual background is modeled with an exponential function while 
partially reconstructed B J/if/X decays require parameterization with a complementary error function. 
All parameters of the residual background PDFs are left free in the fit. 

Figure 1 shows the measured J/il/<t> and invariant-mass spectra with fits overlaid. The fits yield 

Ngo = 6640 ± 100 ^ and = 36290 ± 320 —> Jlt//K*^ signal events. The;y^/dof values 

of the fits are 0.959 for and 0.945 for indicating that both fits describe the data well. 

The additive systematic uncertainties result from the ^ and ^ JfiJ/K*^ signal and back¬ 

ground parameterizations. The contribution from the signal shape parameterization is calculated by vary¬ 
ing the five fixed parameters within ±lcr in a multivariate Gaussian function that takes into account all 
correlations. For non-peaking backgrounds, the exponential function is replaced with a second-order 
polynomial for the and with a second-order polynomial plus an error function for the The differ¬ 
ence in signal yield with respect to the nominal fit is taken as a systematic error. For peaking backgrounds, 
the fixed parameters are varied by ±lcr and the difference with respect to the nominal yield is taken as a 
systematic error. In addition, since S-wave contributions from ^ JIi]jK^K~ and 
decays to ^ and —> JlipK*^ are neglected in the fits, an uncertainty is derived using the AT¬ 
LAS measured contribution of 2.4% [22] for JI and the contribution of 1% for ^ Jltf/K*^ 

derived from the MC simulation. All additive systematic errors are added in quadrature, yielding total 
additive uncertainties of 220 Npo and 650 Ano events. 

The multiplicative systematic uncertainty includes contributions from the relative efficiency and the 
branching fractions of the 0 and K*^ decays. The uncertainty on the relative efficiency is dominated by 
the uncertainty on the (pjK*^ selection (1.2%) which is obtained by varying the fixed fit parameters in the 
(j) and fits by ±lcr and adding all contributions in quadrature. Other uncertainties from the Jlijj selec¬ 
tion (0.2%), reweighting (0.4%), B^^ lifetimes (0.002%) and the contribution due to uncertainties 

in the polarization parameters (0.01%) are negligible. Varying the selection criteria of;y^/dof, Lxy, and 
a gives negligible contributions. Table 1 summarizes the contributions of the additive and multiplicative 
systematic errors. 
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Table 1: Measured and signal yields, the efficiency ratio extracted from simulations, world averages for 
(j) and K*'^ decay branching fractions as well as corresponding systematic uncertainties a on ^ ' 


Observable 

Value 

(T 

Ref. 


6640 ± 100 ± 220 

3.3% 


Ngo 

36290 ± 320 + 650 

1 . 8 % 



0.799 ± 0.001 ± 0.010 

1.3% 


S{(p K+K-) 

0.489 ± 0.005 

1 . 0 % 

[15] 


0.66503 ± 0.00014 

0 . 02 % 

[15] 

Tofal 


4.1% 



From the ratio jN^ after efficiency correction and division by ^ and K*^ decay branching fractions, 
ATLAS measures 


fd ^ 


0.199 ± 0.004(stat) ± 0.008(sys). 


A perturbative QCD prediction [23] yields 


(3) 




0 


■"-0.02V 




where the uncertainties result from the shape parameter cug of the B meson wave function, meson decay 
constants fm, Gegenbauer moments a, in the wave functions of the light vector mesons and the c-quark 
mass. Adding all contributions linearly yields a 7.1% theory error. Using this prediction, the ratio of 
fragmentation fractions is measured to be 

4 - 0.240 ± 0.004(stat) ± O.OlO(sys) ± 0.017(th). (4) 

Jd 

Figure 2 (right) shows the ATLAS fs/fd measurement in comparison with results from LEP [ 6 ], CDF [ 6 , 
7] and LHCb [ 8 , 9]. The ratio fsifd may depend on px and 77 of the B meson, e.g. LHCb observes a px 
but no rj dependence of fsifd [ 8 ]- Figure 2 (left) shows the px dependence of fsifd for ATLAS and that 
of other experiments. To investigate the px and p dependences of fs/fd, the data sample is divided into 
six Px bins in the range 8 GeV < px < 50 GeV and into four p bins for |p| < 2.5 such that the number 
of events in each bin is approximately equal. The fs/fd distributions as a function of px and p have been 
fitted with a uniform (first-order polynomial) distribution yielding fit p-values 0.54 (0.66) and 0.66 (0.49), 
respectively. No significant fs/fd dependence on px and |p| is seen at the present level of accuracy. 

In summary, this Letter reports on the first ATLAS measurement of the ratio of B/] J/ff ^rid B^ 
J/i//K*^ branching fractions multiplied by the ratio of fragmentation fractions fs/fd from which fs/fd 
is determined. The data were produced at the LHC in pp collisions at ^fs = 1 TeV and correspond to 
an integrated luminosity of 2.47 fb“^ This fs/fd measurement, obtained with a new approach, agrees 
with the LHCb [ 8 , 9] results improving the world average considerably. A comparison with the CDF [ 6 , 
7] measurement and the LEP [ 6 ] average confirms fhe universalify of fs/fd- The ATEAS dafa show no 
dependence on px nor on |p| wifhin fhe kinematic range fesfed. 
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Figure 2: (Left) Measurements of fs!fd versus B meson pj for CDF [7], LHCb [8] and ATLAS, where the ATLAS 
data points are plotted at the average px of the events in each bin. The error bars show statistical and systematic 
errors added in quadrature. The LEP ratio, taken from Ref. [6], is plotted at an average px value in Z decays. (Right) 
Measurements of fs! fd (black and blue points with error bars) from LEP [6], CDE [6], LHCb [8, 9] and ATLAS. 
The total experimental error (thin black) is added linearly to the theory error (thick red). The green-shaded region 
shows the HEAG average obtained using the blue points. 
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